On the calculation of the Casimir forces 
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Casimir forces are a manifestation of the change in the zero-point energy of the vacuum caused 
by the insertion of boundaries. We show how the Casimir force can be computed by consideration 
of the vacuum fluctuations that are suppressed by the boundaries, and rederive the scalar Casimir 
effects for a series of geometries. For the planar case a finite universal force is automatically found. 
For curved geometries formally divergent expressions are encountered which we argue are largely 
due to the divergent self-energy of the boundary contributing to the force. This idea is supported 
by computing the effect for a fixed perimeter wedge-arc geometry in two dimensions. 
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Casimir interactions are due to the macroscopic re- 
sponse of the physical vacuum to the introduction of 
boundaries. They were first derived as an attractive 
force between perfectly conductive parallel plates induced 
by the zero-point motion of the electromagnetic field pj 
There is convincing experimental evidence for the reality 
of these forces and a vast body of literature dedicated 
to various aspects of the phenomenon jsj . 

The Casimir interaction £ is the difference between the 
vacuum energy of the system constrained by the bound- 
aries and that of free space. Since boundaries made of 
real materials are transparent to sufficiently high-energy 
modes [1], the high energy spectrum is unaffected by 
the geometry of the system, and only a finite range of 
the spectrum need be considered. However, in the theo- 
retical treatments of this effect the vacuum energies are 
usually calculated from an effective low-energy harmonic 
field theory (such as quantum electrodynamics in the case 
of the electromagnetic Casimir effect), so that they are 
approximated by the sum of zero-point energies of a col- 
lection of simple harmonic oscillators with a spectrum 
uj = c|k| (where c is the speed of light). In this model, 
the dispersion relation holds for arbitrarily large wave 
vectors k; both the "constrained" and "free" vacuum en- 
ergy densities are ultraviolet divergent; and the Casimir 
interaction is the difference between two infinite quanti- 
ties. The theory resolves this problem by a soft-cutoff 
modification of the large-k part of the spectrum that 
leads to a finite vacuum energy. The result is not very 
sensitive to the form of the cutoff, so that the Casimir 
interaction can then be extracted by taking the cutoff 
to infinity at the end of calculation. Other approaches 
to the calculation that make use of analytic continuation 
and dimensional regularization [H techniques give the 
same answer, thus adding to the credibility of the result. 

As already noted, the divergences are more mathemat- 
ical artifacts than physical reality. The important virtue 
of the model is that for many geometries a finite result 
is obtained without introducing a cutoff, leading to a fi- 
nite universal result depending only on h, c, and macro- 
scopic length scales. However, this is not always the case: 



specifically, the divergences occurring for spherical ge- 
ometry in even space dimensions do not cancel [^, 
What this means physically represents an open problem; 
it seems to imply that in two dimensions a conducting 
ring placed in vacuum is unstable. 

The goal of this Letter is two-fold: first, we show how 
the Casimir effect can be computed efficiently by direct 
consideration of the fluctuation modes that are elimi- 
nated by the presence of the boundaries. Second, we in- 
terpret the divergences encountered in the case of curved 
boundaries as being due to divergent self-energy of the 
boundary contributing to the Casimir force. 

These ideas will be illustrated by analyzing Casimir ef- 
fects in a Gaussian field theory with the Euclidian action 
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where T is the temperature and the real scalar u, a func- 
tion of d-dimensional position vector r and imaginary 
time r, is periodic on the Matsubara circle, u(r, 0) = 
u{v, h/T) [8|. The action ([T]) is applicable at energies low 
compared to some scale tiujQ\ corresponding cutoff func- 
tion will be suppressed and invoked only as needed. 

Assume the vacuum is disturbed by the presence of 
sharp boundaries Di and/or Nj of the Dirichlet, wjci = 
0, and/or Neumann, du/dn\i^. — 0, type, respectively, 
where the subscripts i and j label the boundaries and 
du/dn is the normal derivative. The vacuum fluctuations 
eliminated by these constraints are field configurations 
which are solutions to the boundary-value problem for 
the Laplace equation 

(^2972+'^)" = 0' -/.(r,r),— =ffj(r,r) (2) 

where and gj are functions defined on the boundaries 
and playing a role of dynamical variables of our approach. 
The Casimir energy is the negative of the zero-point en- 
ergy of the field configurations satisfying ([2]) as the lat- 
ter do not contribute into the energy of the vacuum dis- 
turbed by the boundaries. Since they are static, a sim- 
plification is achieved by expanding all the dynamical 
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variables of the problem into Fourier series in imaginary 
time domain; for example u(t,t) = ittj(r) exp iwr 
where the Fourier coefficients u^{r) are solutions to the 
boundary-value problem for the Helmholtz equation 

(A-^)w^ = 0, u^\d, = LA^)^ " 5a;j(r) (3) 

Calculations of the zero-point energy of the eliminated 
modes are further simplified if the identity (Vw)^ = 
div{uVu) — uAu is substituted into the action ([T]). Then 
the integral of div{u\/u) over d'^x transforms into a sum 
of surface integrals. The remaining integral over dr van- 
ishes due to the relation Am = —d'^u/c^dr'^ and the con- 
dition of periodicity, M(r,0) = u{r, h/T). As a result we 
find 

rh/T 



1 /"'VJ r 
S'b = - y dT22 J [uVu]i,dsi 

E J [u^Vu-^l^dsi (4) 
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Here [-0] stands for the discontinuity of ip across the 
boundary, and the summation is performed over all the 
boundaries. 

The calculation of the zero-point energy of the elimi- 
nated field configurations takes advantage of the corre- 
spondence between the Feynman path integral for the 
d-dimensional field theory with the action Se, Z = 
J Du{r,T) ex.p{—SE[u]/h), and a partition function of 
a, d + 1-dimensional classical statistical mechanics prob- 
lem with the Hamiltonian Se at a fictitious temperature 
which is equal to Planck's constant The "free en- 
ergy" per unit "length" in the imaginary time direction, 
—h{ln Z)/{H/T) = —TlnZ, gives the zero-point energy 
of the field configurations eliminated by the boundaries. 
Thus the Casimir energy is £ = TlnZ. This completes 
the general formulation of the method which we now il- 
lustrate by analyzing zero-temperature Casimir effects in 
several geometries. We restrict ourselves to the case of 
Dirichlct boundaries. 

Planar geometry. Consider three Dirichlet planes at 
z = 0, z = a, and z = L, where z is one of the axes 
of the d-dimensional rectangular coordinate system and 
< a < L. The outer boundaries are fixed in place so 
that there is no need to look beyond them. We are inter- 
ested in the Casimir pressure exerted on the middle parti- 
tion at z = a. Since the space is uniform relative to trans- 
lations parallel to the boundaries, the field Uui{r) is ex- 
panded into a Fourier series Maj(r) = '^q^'^'ii-^) expiqrj^ 
where r^ is the position vector perpendicular to the 
z axis. Then the boundary- value problem ([3]) for the 
Fourier coefficients Ui^q{z) becomes 



dz2 



-)u^ci = 0, U^q|o,L = 0, U^qla = /wq (5) 



The solution to ([5]) is 

sinh(|K|z) 



sinh(|K|a) ' 



^ z ^ a, K 



_ sinh(|/c|(L-z)) 
smh(|K|(L — ajj 



a < z ^ L 
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Substituting this in Eq. (j4|) we see that only the partition 
at z = a contributes into the action Se with the result 

= ^T.\^\icotH\K\a)+cothi\t,\iL^a)))\Uci\' (7) 

a),q 

where A is the macroscopic {d — l)-dimensional area of 
the boundary. Then up to an a-independent constant the 
Casimir energy per unit area is given by 

f = E 1^ (coth(|«;|a) + coth(|K|(L - a))) 

w,q 
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dujd'^^^q 

— , In (coth(|K|a) -f- coth(|K|(L — a))) 



K dKln(coth(«;a) -|- coth(K(iy — a))) 

(8) 



where in taking the T = and macroscopic lim- 
its we used the rules (f^/T) J duj/2n, J2q ^ 
A j d'^~^q/ {2'kY~^ , respectively. The parameter Kd 
in the third representation is the surface area of a 
rf-dimensional unit sphere, 27r''/2/r(d/2), divided by 
(27r)'*. The Casimir pressure on the boundary, V = 
—d{£/A)/da can be found in closed form 



V = 



dT{^)ad+\] 

(47r)^ 
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where r(a::) and C,{x) are Euler's and Riemann's gamma 
and zeta functions, respectively fioll. In arriving at Q 
we used the gamma function duplication formula [10] and 
the value of the integral 2:''(cotha; — l)dx = 2^''r(d + 

i)C(rf + i) [HI. 

We see that the partition at z = a is attracted to 
the closest outer boundary; specifically, in one dimension 
Eq.® reduces to the well-known result Since the 

outer boundaries impose the Dirichlet boundary condi- 
tions we can imagine joining them together. Then Eq.® 
describes Casimir interaction between two boundaries; 
taking the L — > oo limit we then reproduce the result of 
Ambj0rn and Wolfram 

Circular geometry. Consider a Dirichlet circle of 
radius a in two spatial dimensions. The boundary-value 
problem fS]) for this geometry becomes 



1 d'^ 
p2 dip^ 
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where p and (p are the polar coordinates. Seeking the par- 
ticular solution in the form, u^^ {p, ip) = (p) exp imp, 
where n is an arbitrary integer, we find that the radial 
function Ruj{p) satisfies the equation 



(PR^ ^1 dR^ 



dp^ p dp 
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whose linearly-independent solutions are modified Bessel 
functions In{\Lo\p/c) and Kn{\uj\p/c) [lo| . Thus the gen- 
eral solution to the boundary-value problem (|10p finite 
at p = and decaying as p — > oo is 



IniHa/c) 

Kni\uj\p/C 



n— — oc 



f^n exp irvp, p^a 
fc^nexpinip, p> a (12) 



^ Kn{\uj\a/c 

At the circle p = a both of these reduce to fu;{v) = 
S^-oo f'^n exp inip. Substituting the solution in 
Eq.lgl), using the Wronskian Kn{z)I^{z) - In{z)K^{z) — 
1/z [10], and performing angular integration we find 



l/wri I 



{\u!\a/c)Kn{\uj\a/c) 

The Casimir energy £ and the force J- 
then given by 



(13) 



'd£/da are 
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d/s: In {constIn{tia)Kn{KCi)) (14) 



zdz— (ln(/„(z)if„(z))) (15) 
«^ 



where the const includes all a-independent parameters 
not contributing into the force. The results ([Till and 
(113 are due to Sen [H (see also Refs. @,[ll). 

We have verified that the Casimir effects for spherical 
and cylindrical geometries 15 1 and the cases of the 
Neumann and mixed boundary conditions can be treated 
similarly; unfortunately, we have no space here to enter 
into details of the analysis. 

Divergences and their interpretation. Although 
Eq.® solves the planar version of the problem predict- 
ing finite universal Casimir pressure, its circular (and any 
curvilinear) counterpart (|15p is divergent. Various re- 
searchers dealt with this issue differently. Sen ,13] viewed 
the action ([1]) as an effective low-energy theory which 
must be supplemented by a cutoff function. Introduc- 
ing this into the integrand of Eq. (jl4p removes the diver- 
gence, thus predicting finite non-universal effect. Other 
researchers attempted to remove the divergences by a va- 
riety of techniques [6*] but a divergence was still found in 
the two-dimensional circular case. 
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Figure 1: Wedge of opening angle /3 with superimposed arc 
of radius a in two dimensions. 



In order to understand the difference between planar 
and curved geometries we notice that the Casimir force 
is the change of the energy upon infinitesimal displace- 
ment of the boundary. A geometrically sharp boundary 
possesses a formally divergent energy per unit area. This 
divergent self-energy does not contribute to the Casimir 
force in the planar case as the overall area remains fixed 
as the boundary is displaced. However this is not the 
case for curved boundaries. Indeed a change of the ra- 
dius of a circle implies a change of the perimeter and as 
a result the divergent self-energy will contribute into the 
force. We argue that this is the underlying reason behind 
divergences encountered for curved geometries. In order 
to test this idea we need to look at a curved geometry 
where the area of the boundaries can be kept fixed. 

Wedge-arc geometry in two dimensions. Con- 
sider a wedge of opening angle f3 with superimposed arc 
of radius a. Fig. 1. This is the geometry where both a 
Casimir force and torque occur [16|. Additionally, if the 
variations of (3 and a are constrained by the condition 
of fixed arc length I = Pa, the net perimeter of the in- 
finite wedge edges and the arc are fixed. The Casimir 
interaction here can be inferred from the results for the 
circular geometry since the boundary-value problem we 
need to solve is closely related to (fTO| . The difference 
is that we seek a solution inside the Dirichlet wedge, 
^ (f ^ (3 thus implying Ui^{p,ip) = Ri^{p) sui{'Krnp / (3) , 
11 = 1,2, ... for the particular solution. The radial func- 
tion Ru,{p) satisfies the same Eq.(fTT|) with n being re- 
placed by im/ f3. As a result the solution in question can 
be obtained from Eq. (fT^ by replacing the order n of the 
Bessel functions with nn/P, the angular function exp inip 
with sin(7rn(/j//3), and restricting the summation over n 
from 1 to infinity. The calculation of the Casimir energy 
is similar to that for the circular geometry with the result 



(18) 
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where the const is both a and /3 independent. We then 
employ the uniform asymptotic n ^ 1 expansion of De- 
bye Bill: 



x'^y/^ L,i{nx)Kninx) 
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which can be used to evaluate the energy (I16|) in the 
/3 <C 1 limit. To leading order we find 
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V27rn(l + t2)i/2 



(18) 



Although this expression is formally divergent, it depends 
on a and P only through the fixed arc length I = (3a. 
Thus (fT8|) is a constant which can be subtracted from 
(I16p . As a result we find for presumably universal part 
of the Casimir energy U = £ — 5^^^: 

fc oo -.oo 

he w-^ / - 7rn , ,1 ,i:n , ,7rn 

^E- 1 dtH2j^ii+tr^i^iy)K^iy)) 

(19) 

Using the expansion p7p again we evaluate the energy 
to the next order (fT9|) in the /3 ^ 1 limit with the result 
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which is marginally divergent. This means there is a 
weak dependence on the cutoff frequency wq. The latter 
can be recovered by employing the equality J2n=i "-"^ ~ 
C + lnTV + eAT where C is Euler's constant and cat ^ as 
iV — > (X) The parameter N is estimated by recalling 
that a cutoff function F{cz/u)oa) is suppressed from the 
integrand of Eq.((T6]). In terms of the variable t of Ea. lfTOl) 
this becomes F {imct / uJol) which removes the divergence 
in (|20p by effectively ending the sum of 1/n at n = 
such as TTnet/ujQl « 1 and t ^ 1. Thus A^ « WqI/c. As a 
result we find that with logarithmic accuracy 



hcP'^ ujqI 

U = 7 In , 

2567r/ c 



/3«1, ^»1 (21) 
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This result supports our idea that the divergences are 
due to the boundary self-energy because the dependence 
on the cutoff frequency o^o is logarithmically weak and 
the amplitude of the logarithmic dependence is universal. 
Sen's result for the circle JJi] also has a term logarithmi- 
cally dependent on the cutoff frequency ujq. 

Our result (|2T|) means that for fixed arc length I the 
zero-point motion induces a widening torque —dU/dl3 
which is not unexpected on the physical grounds. The ac- 
curacy of the Debye expansion (|17p , implies that Eq. (I^T]) 
approximately captures the whole ^ /3 ^ 27r range. 



To summarize, we have demonstrated how Casimir ef- 
fects caused by sharp boundaries can be efficiently com- 
puted by focusing on the quantum fluctuations elimi- 
nated by these boundaries. The applicability of this 
method is not limited to the scalar field theory ([1]), 
Dirichlet boundaries, zero-temperature limit or to the ge- 
ometries we have considered. 

Additionally, we argued that the formally divergent 
Casimir forces encountered in the presence of curved 
boundaries are due to divergent self-energy contributions, 
and supported this idea by an explicit calculation. More 
work is needed to further test this hypothesis. Our anal- 
ysis also supports Sen's viewpoint [13[ that the Casimir 
effect in a Dirichlet ring in two dimensions is finite and 
non-universal with the cutoff frequency wq supplied by 
the properties of the material the boundary is made of. 

This work was supported by the Thomas F. Jeffress 
and Kate Miller Jeffress Memorial Trust. 
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